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abstract 

A  technique  is  described  for  using  recursive  digital  filters  in  the 
moving  target  indicator  of  phased  array  radars.  This  technique,  called 
initialization,  is  effective  in  reducing  the  severe  transient  effects 
normally  exhibited  by  these  filters  for  large  clutter  input  signals. 
Consequently,  the  superior  frequency  characteristic  of  recursive  filters 
can  be  realized.  The  problem  is  addressed  from  the  state  variable  point 
of  view  and  a  general  initialization  expression  is  derived.  This  pro¬ 
cedure  is  illustrated  for  example  filters ,  and  typical  moving  target 
indicator  output  sequences  are  shown.  The  concept  of  "dynamic"  frequency 
response  is  introduced  and  used  to  demonstrate  the  effectiveness  of 
initialization  for  moving  target  indicator  input  sequences  of  varying 
length . 
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1.  Introduction 

Many  radars  presently  employ  a  digital  implementation  of  the 
delay-line  canceler  for  the  detection  of  moving  targets  in  a  background 
of  clutter.  This  moving  target  indicator  (MTI),  which  is  actually  a 
nonrecursive  digital  filter,  may  be  reasonably  effective  against  low 
frequency  components  of  clutter.  The  standard  cancelers,  however, 
exhibit  rather  poor  pass-band  frequency  characteristics,  and  a  number 
of  techniques  have  been  developed  to  improve  their  response  by  using 
different  feed-forward  multiplier  coefficients  |1).  These  nonrecursively 
implemented  filters,  which  possess  finite-duration  impulse  responses, 
are  somewhat  limited  in  the  amount  of  frequency  response  shaping  that 
can  be  obtained.  Recursively  implemented  filters,  on  the  other  hand, 
offer  superior  frequency  characteristics  over  nonrecursive  filters  of 
comparable  hardware  [2).  Due  to  the  feedback  nature  of  these  filters, 
however,  the  transient  response  due  to  an  impulse  input  is  much  more 
severe.  In  fact  these  recursively  implemented  filters  are  frequently 
called  infinite-duration  impulse  response  filters. 


2.  Transient  Problem  for  Recursive  MTI  Filters 

Because  of  the  relatively  poor  transient  performance  of  recur¬ 
sive  filters,  these  filters  have  seldom  been  employed  for  MTI  purposes. 
T.e  transient  response  problem  is  particularly  severe  in  a  strong 
ground  clutter  environment.  For  this  case  the  clutter  returns  may  be 
orders  of  magnitude  larger  than  those  of  a  target  and  effectively  appear 
as  a  large  step  input  to  the  filter.  This  can  produce  such  excessive 
ringing  in  the  filter  output  that  the  target  will  be  masked  until  the 
transient  response  has  settled.  Furthermore,  the  oscillatory  charac¬ 
teristic  of  the  MTI  increases  in  magnitude  as  the  bandwidth  of  the 
filter  increases,  i.e.,  as  the  clutter  notch  at  zero  doppler  frequency 
is  narrowed.  This  effect  is  graphically  depicted  In  Figure  1  which 
shows  the  unit  step  response  for  three  four-pole  Butterworth  filters. 

As  the  3-dB  cutoff  frequency  is  lowered,  the  transient  response  of  these 
filter®  becomes  increasingly  severe.  Sluvc  uiu&l  ludars  are  limited  in 
the  time  allocated  to  each  sector  of  their  scan,  the  long  settling 
times  required  reduce  the  effectiveness  of  the  recursive  implementation. 

For  conventional  scanning  radars  the  modulation  due  to  the  scanning 
antenna  pattern  reduces  the  severity  <,f  the  transient  response  from 
large  sharply  defined  clutter  sources.  Nevertheless  the  transient 
problem  may  still  be  significant,  and  a  trade-off  between  frequency 
response  characteristics  and  acceptable  transient  behavior  is  usually 
required.  While  this  has  resulted  in  some  satisfactory  recursive  MTI 
designs,  the  nonrecursive  implementation  has  generally  been  employed. 
However,  for  radars  that  use  step  scanning  such  as  those  with  phased 
array  antennas,  a  new  technique  [3}  has  been  suggested  to  minimize  the 
poor  transient  performance  of  recursive  MTI  filters. 
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3.  Initialization  Technique 

Thts  technique,  which  is  called  initialization,  takes  advantage 
of  the  discrete  nature  of  phased  array  scanning  to  reduce  transient 
effects.  In  contrast  with  conventional  scanning  systems  which  process 
returns  continuously,  phased  array  systems  reinitiate  the  processing 
each  time  the  beam  is  stepped  to  a  new  position.  Hem  e,  the  initial 
return  from  each  new  beam  position  can  be  uaed  to  apply  initial  condi¬ 
tions  to  the  MTI  for  processing  the  remaining  returns  from  that  position. 
Since  ground  clutter  usually  exhibits  relatively  low  frequency  compo¬ 
nents,  the  clutter  returns  can  be  approximated  by  a  step  input  equal  in 
magnitude  to  the  first  return  for  each  beam  position.  The  steady-state 
values  that  would  normally  appear  in  the  filter  memory  i'lements  after  an 
infinitely  long  sequence  of  these  inputs  can  therefore  he  immediately 
calculated  and  loaded  into  the  filter  to  suppress  the  transient  response. 

The  initialization  technique  can  be  derived  by  considering  the 
general  state  variable  block  diagram  of  Figure  2.  This  linear,  time- 
invariant  single  input-single  output  system  represents  the  digital  MTI 
with  the  state  variables  x(n)  defined  as  the  outputs  of  the  MTI  filter 
registers .* 


Figure  2.  General  State  Variable  Representation  of 
Digital  MTI 

The  corresponding  state  variable  equations  are  [4] 


x(n  +  1)  *  Ax(n)  +  Bu(n) 

(1) 

y(n)  -  Cx(n)  +  Du(n) 

(2) 

*The  abbreviated  notation  jt(n)  A  Jt(nT)  is  used  in  this  work. 
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and  the  solution  to  the  state  equation  is 


n-  i 

x(n)  =  Anx(0)  +  ^  An-1_mBu(m)  .  (3) 

m=0 


For  the  initialization  technique  under  consideration,  the  basic  pro¬ 
cedure  is  to  immediately  force  the  registers  of  the  filter  to  their 
steady-state  values  by  applying  initial  conditions  to  the  filter. 
Therefore,  the  required  initial  conditions,  x(0),  must  be  such  that  for 
a  step  input  the  states  remain  constant,  specifically, 


x(n  +  1)  =  x(n)  =  x(0) ,  for  all  n  >  0  .  (4) 


It  then  follows  from  (3)  that 


x(n  +  1)  -  x(n)  *  (A1*1  -  An)  x(0)  +  A^BuCO) 

n- 1  n- 1 

+  ^  An_1"mBu(m+  1)  -  ^  An"l“mBu(m)  .  (5) 

n>“0  n*>0 

Substituting  x(n  -f  1)  -  x(n)  “  0  from  (4)  and  u(m  +  1)  -  u(m)  =  0  from 
the  step  input  assumption  yields 

0  -  An(A  -  I)x(O)  +  AnBu(0)  .  (6) 


Consequently,  the  initialization  expression  can  be  written  as 


x(0)  -  (I  -  A)-1Bu(0) 


(7) 


by  excluding  the  trivial  case  ^  »  0  and  assuming  that  the  inverse  in 
(7)  exists.  This  expression  gives  the  values,  x(0)»  that  must  be  stored 
in  the  MTI  memory  registers  once  the  initial  return  from  each  new  beam 
position,  u(0),  is  known. 

Although  the  use  of  (7)  will  result  in  the  suppression  of  the 
transient  response  to  a  step  input,  it  does  not  ensure  that  the  MTI 
output  will  be  at  a  zero  level  in  its  steady-state  condition.  This 
suggests  that  a  restriction  should  be  placed  on  the  two  coupling 
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matrices  C  and  D  of  the  output  equation  (2)  to  ensure  that  the  steady- 
state  output  to  a  step  input  is  zero.  The  restriction  can  be  obtained 
by  substitution  of  (3)  into  (2)  giving  an  expression  for  the  MTI  output 
as  a  function  ot  the  input  and  initial  states,  i.e., 


y (n)  =  C 


Anx(0) 


n-  1 

s* 

n*»0 


p-  L- 


Bu  (o) 


+  Du(n) 


(8) 


Now,  by  using  the  initialization  expression  of  (7),  the  output  can  be 
written  as 


n-1 

y (n)  =  C  An(I  -  A)_1Bu(0)  +  C  J  An'1-mBu(m)  +  Du(n)  .  (9) 

iw=0 

For  a  step  input  it  is  desired  that  the  output,  be  zero  for  all  n  and  in 
particular  for  n  =  0.  Consequently,  (9)  reduces  to 

0  =  C  A°(I  -  A) ”  1JBu(0)  +  Du(0)  (10) 

and  it  follows  that  the  desired  relationship  is 

D  -  -C(I  -  A)-lB  .  (11) 


4.  Example  MTI  Filters 

The  application  of  (7)  and  (LI)  can  be  illustrated  by  consider¬ 
ing  the  MTT  ftlter  of  Figure  3(a).  By  inn  pec  Lion  of  this  block  diagram, 
the  state  variable  equations  are 


and 


xL(n  +  L ) 

"0  1  * 

(n)] 

"o' 

- 

+ 

x2(n  +  i) 

_b2  \ 

_x2(n)J 

1 

u(n) 


(12) 


y(n)  -  [l  +  b2  j  bL 


’] 


"*l  (n)* 
-X2(n)- 


+  fl)u(. 


(13) 
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It  then  follows  from  (7)  that  the  initialization  values  to  be  used  are 


[*L<0)" 

L 

"o' 

ufm  u(0) 

T 

Lv°>. 

“  1  '  bl  "  b2 

b2  l. 

1 

(0)  "  1  -  bx  -  b. 

l 

(14) 


and  from  (11)  the  required  D  is 


D  = 


1  -f  b. 


{ 1 i  (15) 


which  is  satisfied  for  Figure  3(a).  Consequently,  the  first  return  from 
each  beam  position  is  multiplied  by  (1/(1  -  b^  -  b^)  ],  this  product  is 
Loaded  in  the  filter  registers,  and  the  remaining  returns  from  that  beam 
position  are  processed  in  the  normal  manner.  However,  (14)  illustrates 
a  problem  with  the  conventional  implementation  scheme  of  Figure  3(a). 

As  the  filter  cutoff  frequency  approaches  zero,  the  denominator  of  the 
me f f i r f on f  term  in  (14)  alee  approaches  zero.  This  i’Csults  In  u  vciy 
large  magnification  of  the  input  signal  and  a  corresponding  increase 
in  the  number  of  bits  required  to  represent  the  register  values.  An 
Improved  implementation  [3]  having  the  same  steady-state  frequency 
characteristic  as  the  conventional  implementation  but  different  state 
variable  equations  is  shown  in  Figure  3(fc).  Since  the  state  equations 
a  re 


(n  +  1) 

*0  l  ' 

(n) 

-1 

ta 

+ 

x2(n  +  1) 

_b2  bL. 

x?(n) 

Wm  ^  M 

(16) 
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and 


(17) 


the  initialization  values  are 


x2(°) 


f-  - 

-1 

r-  - 

0 

u(0)  = 

}  -  bl. 

u(0)_ 

(18) 


and  the  required  D  is 


ri  -  b  ii 

'  -1  ' 

-  -l  fl  *  ll 

L 

1  -  bj  -  b  [L  :  LJ 

b.  1 

1  -  b, 

L  2  J 

L  lJ 

[1]  •  (19) 


Expression  (18)  indicates  a  particularly  convenient  form  of  initial¬ 
ization  in  which  the  first  memory  register  of  the  filter  is  simply  loaded 
with  the  initial  radar  return  while  the  second  register  is  zeroed.  The 
registers  in  any  additional  cascaded  sections  following  the  first  are 
also  zeroed.  Moreover,  the  dynamic  range  requirements  for  the  registers 
of  the  modified  implementation  are  reduced  in  comparison  with  those  of 
the  conventional  implementation  since  "he  multiplier  coefficient  in  (18) 
if.  now  unity.  Both  of  these  implementations  satisfy  the  relationship  of 
(11)  which  ensures  that  the  steady-state  output  will  be  zero  for  a  step 
input.  Any  nc~«l  MTI  filt^t  wiui  a  clutter  notch  ot  zero  gain  at  zero 
frequency  will  satisfy  this  relationship. 


5.  Four-Pole  Filter 

The  filters  of  the  preceding  section  may  be  cascaded  as  shov*n 
in  Figure  4  to  give  a  very  useful  filter  configuration.  As  previously 
mentioned,  the  initialization  of  this  filter  consists  of  loading  u(0) 
into  the  first  register,  x^(n)  in  Figure  4,  and  then  zeroing  Che  remain¬ 
ing  registers.  This  can  readily  be  verified  from  the  state  and  output 
equations  for  this  filter,  (20)  and  (21). 
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Figure  4.  Four- Pole  Recursive  Filter 


*1('n  +1) 

0  10  0 

xL(n) 

0 

x2(n  +  1) 

b4  b3  1  -1 

x2(n) 

1 

= 

+ 

x3(n  +  1) 

0  0  0  1 

x3(n) 

-1 

x4(n  +  1) 

0  0  b2  bL 

Vn) 

1  -  bj 

- 

•  - 

- 

_  _ 

u(n) 


(20) 


and 


y(n)  =  [l  +  b4  •  aL  +  b3  •  1  •  '  L] 


xL(n) 


*2(n) 


x3(n) 


x4(n) 


+  El]u(n) 


(21) 


Using  (7)  and  A  and  B  from  (20),  the  initialization  values  are 
seen  to  be 


xL(0) 

0 

x2  (0) 

0 

x,(0) 

m 

0 

x4(°) 

u(0) 

—  ota 

which  verifies  the  preceding  comments. 


(22) 
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Another  important  feature  of  this  filter  is  that  (LL)  is  satisfied 
regardless  of  the  value  of  a^.  Therefore,  the  first  filter  section 

assures  a  zero  dc  gain  of  the  filter,  and  the  of  the  second  section 

can  be  used  to  provide  a  notch  in  the  frequency  response  at  some  fre¬ 
quency  other  than  zero. 


6.  Transient  Response  Comparisons 

It  has  been  assumed  for  the  purposes  of  initialization  that 
the  first  return  is  due  only  to  clutter,  but  obviously  this  is  not 
always  the  case.  If  a  target  signal  is  present  along  with  the  clutter, 
there  may  be  an  error  in  the  initialization  value  which  will  produce  a 
transient,  effect  in  the  filter  output.  A  similar  situation  exists  when 
there  is  a  target  signal  alone.  Generally,  these  effects  are  small  when 
compared  with  the  ringing  caused  by  ciutter  in  a  recursive  1  in  with  no 
initialization. 

The  improvement  in  typical  MTI  transient  responses  is  illustrated 
in  Figure  5  for  systems  with  and  without  initialization.  These  response 
are  for  a  our-pole  Bucterworth  filter  having  a  3-dB  cutoff  frequency  of 
4  percent  of  the  radar  pulse- repetition  frequency  (PRF).  The  unit  step 
response  of  this  filter  was  given  in  Figure  L (b) .  The  power  spectral 
density  of  the  simulated  ground  clutter  return  consisted  of  a  zero- 
frequency  component  and  a  Gaussian  fluctuating  component  with  an  rms 
f. euaency  spread  equal  to  0.3  percent  of  the  FRF  and  ratio  of  stationary 
to- flue tua ting  clutter  power  of  0.8  [6,  7] .  The  plots  of  Figure  5(a) 
are  the  MTI  responses  for  a  clutter  return  only  while  Figure  5(b)  shows 
tin  filter  output  when  the  input  signal  consists  of  clutter  plus  a 
target.  The  target  return  had  a  doppler  frequency  equal  to  L0  percent 
of  the  PRF  and  the  c lu l ter- to- targe t  power  ratio  was  25  dK  .  It  is 
apparent  from  the  responses  without  initialization  that  the  transient 
effect  marks  the  desired  target  signal  for  an  appreciable  period  of 
time.  With  initialization,  on  cFe  other  hand,  this  effect  is  signifi¬ 
cantly  reduced  and  the  target  can  be  immediately  detected. 


7.  Dynamic  Frequency  Response 

The  effectiveness  of  the  initialization  technique  for  other 
doppler  frequencies  was  investigated  by  simulating  a  radar  quadrature 
processor  and  obtaining  a  type  of  frequency  response  which  incorporates 
the  transient  effects.  The  processor  input  signal  consisted  of  the  in- 
phase  (I)  and  quadrature  (Q)  components  of  a  complex  signal  cxp(j2nf^T) . 

For  each  doppler  frequency,  ,  chosen  in  the  Lange  0  to  DRF/2,  input 

sequences  of  N  pulses  were  generated  to  represent  these  1  and  Q  signals. 
These  sequences  were  processed  by  the  respective  MTl's  and  then 
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recombined  by  taking  the  square  root  of  the  sum  of  tne  squares  of  the 
Mil  outputs.  The  N  recombined  values  were  then  integrated  to  give  a 
processor  output  residue  3(f,).  This  process  is  illustrated  in 
Figure  6. 


Figure  6.  Radar  Quadrature  Processor 


Figures  7  through  9  show  plots  of  the  normalized  S(fd)  values  both 

with  and  without  initialization  for  pulse  sequences  of  varying  lengths. 
These  figures,  called  "dynamic"  frequency  responses  [3],  are  for  the 
three  Butterworth  filters  whose  step  responses  are  given  in  Figure  1. 
These  curves  clearly  indicate  the  improved  performance  offered  by  the 
initialization  technique.  It  is  also  apparent  from  these  curves  that 
the  "dynamic"  frequency  response  characteristic  approaches  the  steady- 
state  MTI  response  as  the  number  of  pulses  processed  increases. 


An  initialization  technique  has  been  described  which  provides 
significant  improvement  in  the  MTI  performance  of  phased  array  radars. 
Moreover,  this  improved  performance  is  obtained  with  virtually  no 
increase  in  hardware  over  a  conventional  recursively  implemented  filter. 
It  has  been  demonstrated  that  the  improvement  is  a  function  of  the  number 
of  pulses  per  beam  position  and  is  degraded  as  the  number  of  pulses  is 
decreased.  However,  for  as  few  as  six  pulses  per  position,  the  initial¬ 
ized  recursive  filter  appears  to  have  an  adequate  characteristic.  It 
should  be  emphasized  that  these  curves  were  obtained  using  the  coeffi¬ 
cients  for  a  standard  Butterworth  filter  which  is  optimum  only  for 
steady-etate  performance.  A  further  improvement  should  be  obtained  if 
MTI  multiplier  coefficients  are  used  that  optimize  the  processor  dynamic 
response  for  a  given  number  of  pulses. 


1.1 


WITHOUT  INITIALIZATION 


V«*F 


Figure  7.  Dynamic  Frequency  Responses  £or  Filter  Cutoff 
Frequency  of  0.02  PRF 
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WITHOUT  INITIALIZATION 


Figure  0.  Dynamic  Frequency  Responses  for  Filter  Cutoff 
Frequency  of  0.04  PRF 
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